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ABSTRACT 

We present a simple, semi-analytical model to compute the mass functions of dark 
matter subhaloes. The masses of subhaloes at their time of accretion are obtained 
from a standard merger tree. During the subsequent evolution, the subhaloes experi- 
ence mass loss due to the combined effect of dynamical friction, tidal stripping, and 
tidal heating. Rather than integrating these effects along individual subhalo orbits, 
we consider the average mass loss rate, where the average is taken over all possible 
orbital configurations. Under the ansatz that the average distribution of orbits is in- 
dependent of parent halo mass, this allows us to write the average mass loss rate as a 
simple function that depends only on redshift and on the instantaneous mass ratio of 
subhalo and parent halo. After calibrating this model by matching the subhalo mass 
function (SHMF) of cluster-sized dark matter haloes obtained from high-resolution, 
numerical simulations, we investigate the predicted mass and redshift dependence of 
the SHMF. We find that, contrary to previous claims, the subhalo mass function is 
not universal. Instead, both the slope and the normalization depend on the ratio of 
the parent halo mass, M, and the characteristic non-linear mass M*. This simply re- 
flects a halo formation time dependence; more massive parent haloes form later, thus 
allowing less time for mass loss to operate. We predict that galaxy-sized haloes, with 
a present-day mass of M ~ lO 12 /^ 1 M have an average mass fraction of dark matter 
subhaloes that is a factor three lower than for massive clusters with M ~ 10 15 /i -1 M . 
We also analyze the halo-to-halo scatter in SHMFs, and show that the subhalo mass 
fraction of individual haloes depends most strongly on their accretion history in the 
last ~ 1 Gyr. Finally we provide a simple fitting function for the average SHMF of 
a parent halo of any mass at any redshift and for any cosmology, and briefly discuss 
several implications of our findings. 

Key words: galaxies: halos — cosmology: theory — dark matter — methods: sta- 
tistical 



1 INTRODUCTION 

During the hierarchical assembly of dark matter haloes, the 
inner regions of early virialized objects often survive accre- 
tion onto a larger system, thus giving rise to a population 
of subhaloes. This substructure evolves as it is subjected to 
the forces that try to dissolve it: dynamical friction, tidal 
forces, and impulsive collisions. Depending on their orbits 
and their masses, these subhaloes therefore either merge, 
are disrupted or survive to the present day. 

To fully describe, in a statistical sense, the non-linear 
distribution of mass in the Universe, it is essential that 
halo substructure is taken into account. After all, galaxies 
are thought to reside at the centers of dark matter haloes, 
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which includes dark matter subhaloes. When building a co- 
herent picture of galaxy formation or of galaxy clustering, 
it is therefore of paramount importance that halo substruc- 
ture is taken into account. In particular, we need an ac- 
curate description of the conditional subhalo mass func- 
tion, n(m\M)dm, which gives the number of subhaloes with 
masses in the range m ± dm/2 that reside in a parent halo 
of mass M. Combined with the (parent) halo mass function, 
n(M)dM, this then provides a complete, statistical descrip- 
tion of the abundance of dark matter haloes down to the 
level of subhaloes. In addition, a comparison of n(m\M)dm 
with the conditional luminosity function, <&(L\M)dL (Yang, 
Mo & van den Bosch 2003; van den Bosch, Yang & Mo 2003) 
will yield important insights into galaxy formation and allow 
for a detailed study of galaxy bias. 

Only since a couple of years numerical simulations of 
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structure formation have reached the mass and force resolu- 
tion to allow for a detailed study of dark matter substructure 
(e.g., Tormen 1997; Tormen, Diaferio & Syer 1998; Moore 
et al. 1998, 1999; Klypin et al. 1999a,b; Ghigna et al. 1998, 
2000; Stoehr et al. 2002; De Lucia et al. 2004; Diemand, 
Moore & Stadel 2004; Gill et al. 2004a,b; Gao et al. 2004; 
Reed et al. 2004; Kravtsov et al. 2004). Most of these stud- 
ies have found that in terms of their substructure properties, 
dark matter haloes are homologous; the internal structure 
of a galaxy-sized halo looks just like a rescaled version of 
that of a rich cluster. This would imply that the subhalo 
mass function is independent of parent halo mass. However, 
most of these results are based on small numbers of individ- 
ual haloes, while halo-to-halo variations are expected to be 
fairly large. Combined with uncertainties due to numerical 
resolution and the identification of dark matter subhaloes, 
this means that the statistical significance of these results 
is still unclear. For example, Gao et al. (2004), analyzing 
a relatively large sample of dark matter haloes extracted 
from a large, high-resolution simulation, find that the nor- 
malization of the SHMF depends on parent halo mass. In 
particular, they claim that more massive haloes contain a 
larger mass fraction in subhaloes. Similar results have been 
obtained by Diemand et al. (2004; their Fig. 7) and Kang et 
al. (2004; their Fig. 2). Such a parent halo mass dependence 
might be expected from the fact that more massive haloes 
form later (e.g., van den Bosch 2002), thus leaving less time 
for mass loss to operate. 

Recently, there have also been a number of analyti- 
cal studies of dark matter subhaloes based on the extended 
Press-Schechter (EPS) formalism (Bond et al. 1991; Bower 
1991; Lacey & Cole 1993). Although the EPS formalism 
only yields information regarding parent haloes, it is a log- 
ical next step to simply associate the progenitor haloes of 
a given parent halo (whose properties can be computed us- 
ing EPS) with its present day subhaloes (Fujita et al. 2002; 
Sheth 2003). This, however, ignores the fact that subhaloes 
experience significant amounts of mass loss. A more realis- 
tic approach, therefore, needs to combine this EPS based 
formalism with an analytical description of the (mass) evo- 
lution of dark matter subhaloes. 

Oguri & Lee (2004), following up on a previous study 
by Lee (2004) , presented a semi-analytical model to compute 
the SHMF from EPS, taking detailed account of dynamical 
friction and tidal stripping. They predict that the SHMF is 
virtually independent of parent halo mass. However, an ob- 
vious downside of their approach is that they use the present 
day mass of the parent halo when computing the impact of 
dynamical friction. In reality, the parent halo mass evolves, 
which should have been taken into account (see e.g., Taffoni 
et al. 2003; Zhao 2004). Therefore, it seems likely that Oguri 
& Lee have underestimated the impact of dynamical friction, 
and, since the mass accretion history depends on mass, may 
not have correctly predicted the mass-dependence of the 
SHMF. Zentner & Bullock (2003; hereafter ZB03) and Tay- 
lor & Babul (2004; hereafter TB04) improve on this by in- 
tegrating orbits in the changing potential of the parent halo 
(whose mass accretion history is computed using detailed 
merger trees). By including detailed analytical descriptions 
of dynamical friction and tidal heating and stripping, these 
authors provide detailed, realistic models for the evolution 
of dark matter substructure. Unfortunately, TB04 refrain 



from a discussion of predictions regarding the SHMF, while 
ZB03 only investigate the cosmology-dependence, but not 
the parent halo mass dependence. 

In this paper we follow a similar approach, except that 
we treat the actual mass loss of subhaloes in a very sim- 
ple manner. We only consider the average mass loss rate, 
where the average is considered to be taken over all orbital 
configurations. This means that we do not have to integrate 
individual orbits, and allows us to write the mass loss rate 
as a function of the mass ratio of subhalo to parent halo 
only. Rather than attempting to obtain an estimate of this 
average mass loss rate from first principles, we simply adopt 
a functional form, and adjust the free parameters to match 
the SHMF obtained from numerical simulations. As in ZB03 
and TB04, we take detailed account of the fact that while 
the subhalo looses mass, the parent halo gains mass due to 
its hierarchical growth. After calibrating the model against 
numerical simulations, we use it to investigate the parent 
halo mass and redshift dependence of the SHMF, as well 
as the halo-to-halo scatter. We show that our simple model 
predicts that (i) more massive haloes have a larger mass 
fraction of substructure, (ii) the halo-to-halo scatter is large, 
(iii) the abundance of subhaloes per unit parent halo mass 
is independent of parent mass, and (iv) the subhalo mass 
fraction is larger at higher redshifts. These findings are in 
excellent agreement with the numerical simulations of Gao 
et al. (2004). The main advantage of our model over either 
numerical simulations or the more detailed models of ZB03 
and TB04 is its shear simplicity and computational speed 
that allows a detailed investigation of the dependence of the 
SHMF on cosmology, parent halo mass, and redshift. In ad- 
dition, it provides a simple description of the average mass 
loss rate of dark matter subhaloes, which may be useful, for 
example, to describe the evolution of the mass-to-light ratio 
of satellite galaxies. 

This paper is organized as follows. In Section 2 we give 
a brief overview of the SHMFs obtained from numerical 
simulations. Section 3 describes our method for computing 
SHMFs based on a combination of EPS and a simple model 
for the average mass loss rate of subhaloes. Sections 4, 5, 
and 6 discuss the mass-dependence, the halo-to-halo vari- 
ance, and the redshift dependence of the SHMF, respec- 
tively. In Section 7 we provide a simple analytical fitting 
function for the average SHMF of a halo of given mass and 
redshift. We summarize our results in Section 8. 

Throughout we use m and M to denote the masses of 
the subhalo and the parent halo, respectively. Here the par- 
ent halo mass is defined as the total mass (including that 
of all subhaloes) within a sphere of density 200 times the 
critical density at redshift zero. For brevity we use tp to 
indicate the mass ratio m/M, and we consider it under- 
stood that m, M, and t[> all depend on time, without having 
to write this time-dependence explicitly. A subscript zero 
is used to indicate the present day value (i.e., at redshift 
zero). Unless specifically stated otherwise, we adopt a flat 
ACDM 'concordance' cosmology with 57 m = 0.3, £7a = 0.7, 
h — H a /(W0 kms -1 Mpc -1 ) = 0.7 and with initial density 
fluctuations described by a scale-invariant power spectrum 
with normalization erg = 0.9. 
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2 SUBHALO MASS FUNCTIONS FROM 
NUMERICAL SIMULATIONS 

As discussed in the introduction, numerous studies have 
determined SHMFs from high-resolution numerical simu- 
lations. In Fig. 1 we compare the subhalo mass functions 
from three independent studies (all based on the same 
ACDM concordance cosmology). The solid dots with er- 
rorbars (Poissonian) indicate the average SHMF obtained 
using the simulations described in Tormen, Moscardini & 
Yoshida (2004) from 17 clusters with masses in the range 

3 x 10 14 /i _1 M Q < M < 1.7 x lO 15 /^ 1 M Q . These high- 
resolution simulations were obtained using the technique 
of re-simulating, at much higher resolution, a region of in- 
terest selected from a large cosmological volume. De Lu- 
cia et al. (2004) studied a similar set of 11 high resolution 
re-simulations of galaxy clusters with masses in the range 
7 x lO 13 /^ 1 M Q < M < 1.8 x 10 15 /i _1 M . The dashed line 
in Fig. 1 corresponds to dn/dln(m/M) = 0.016(m/M)~°' 94 , 
which is the power-law relation that best fits the aver- 
age SHMF of this set (obtained by fitting their results by 
eye) J. Finally, the solid line indicates the power-law SHMF, 
dn/dln(m/M) = 0.017(m/M)~°' 91 , that best fits the results 
of Gao et al. (2004) , obtained by fitting-by-eye the average 
SHMF of their 15 haloes with 3 x lO 14 ^ 1 M Q < M < 
W 15 ^ 1 Mq. 

All three SHMFs are in good agreement with each other, 
both in terms of the slope at small m/M and in terms of 
the normalization. In the range —3.5 < log(m/M) < —2.5, 
where the results are most accurate, all three SHMFs agree 
with each other at better than 20 percent. Given the differ- 
ent force and/or mass resolutions of the various simulations, 
and the different techniques used to identify subhaloes, this 
level of agreement is in fact better than what one might 
naively expect. Especially since relatively small samples of 
haloes have been used, which, if halo-to-halo scatter is large, 
may cause significant scatter in these averages. In Section 3.2 
below we will use these SHMFs to calibrate our model for 
the subhalo mass loss rate. 



3 SUBHALO MASS FUNCTIONS FROM 
MERGER TREES 

The aim of this paper is to develop an algorithm that allows 
a fast and reliable computation of subhalo mass functions. 
As discussed in Section 1 two ingredients are essential: a 
method to compute progenitor haloes, and a proper treat- 
ment of the mass evolution of subhaloes. For the former, we 
use a standard merger tree, which we construct using the N- 
branch method with accretion (Somerville & Kolatt 1999; 
hereafter SK99). We adopt the same time-stepping as in 
SK99, and introduce a lower-mass cut-off of mu m = W~ 4 Mo, 
which reflects the effective mass resolution of our merger 
trees. This cut-off is required since the number of progeni- 
tor haloes diverges as the mass goes to zero. Following SK99, 
any mass contained in haloes below the resolution limit is 

t The normalization of the mass functions shown in panel (f) 
of Fig. 1 of De Lucia et al. (2004) is incorrect and needs to 
be translated in the ^-direction by +0.6, (De Lucia, private 
communication) 
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Figure 1. Comparison of SHMFs of parent haloes in the mass 
range 10 14 /i _1 M Q < M < 10 15 h _1 M Q , as obtained by different 
authors using different high resolution numerical simulations. 

accounted for by referring to it as 'accreted' mass (for which 
the prior mass accretion history is not followed back in time). 

A proper treatment of the mass evolution of the sub- 
haloes is more complicated. A subhalo moving on a fixed 
orbit in a static halo experiences mass loss due to tidal 
stripping and heating. If the orbit were to remain fixed, and 
in the absence of tidal heating, the mass loss rate would 
rapidly decline with time, as all mass beyond the tidal ra- 
dius would be stripped after at most a few orbital periods. 
In reality, however, tidal heating continues to 'push' stars 
beyond the tidal radius, where they can be stripped, while 
the orbit evolves due to dynamical friction which causes the 
tidal radius to shrink. Both effects significantly prolong the 
duration and increase the amount of mass loss, which may 
eventually lead to the complete disruption of the subhalo. 
For detailed numerical simulations of the mass loss of dark 
matter subhaloes see Hayashi et al. (2003) and Kazantzidis 
et al. (2004). 

To properly account for the above mentioned effects, 
which depend strongly on the orbital eccentricity (e.g., 
Colpi, Mayer & Governato 1999; Gnedin, Hernquist & Os- 
triker 1999; Taylor & Babul 2001, 2004; Taffoni et al. 2003), 
requires a detailed integration over all individual subhalo 
orbits. This is complicated by the fact that the mass of the 
parent halo evolves with time. If the mass growth rate is 
sufficiently slow, the evolution may be considered an adi- 
abatic process, thus allowing the orbits of subhaloes to be 
integrated analytically despite the non-static nature of the 
background potential (this principle is exploited in the mod- 
els of ZB03 and TB04). In reality, however, haloes grow hi- 
erarchically through (major) mergers, making the actual or- 
bital evolution highly non-linear. 

In order to sidestep these difficulties, we consider the 
average mass loss rate of dark matter subhaloes, where the 
average is to be taken over the entire distribution of orbital 
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Figure 2. Illustration of the processes of merging and evolution. 
Each time-step, such as the step between t = t\ and t = t<i 
shown here, the parent halo mass grows by merging (assumed 
instantaneous at t = ti), while the subhalo masses evolve due to 
mass loss. See text for a detailed description. 

configurations. This removes the requirement to actually in- 
tegrate individual orbits, resulting in a subhalo mass loss 

rate that depends only on the mass ratio ip = m /M*. This 
is true as long as the average distribution of orbital eccen- 
tricities of subhaloes does not depend on parent halo mass. 
The eccentricity distribution, P(e), depends on the orbital 
anisotropy and on the density distribution of the parent halo 
(van den Bosch et al. 1999). Although more massive haloes 
are less concentrated on average (Navarro, Frcnk & White 
1997), which could give rise to a mass dependence of P(e), 
van den Bosch et al. (1999) have shown that P(e) depends 

t Throughout this paper we ignore the fact that the average halo 
concentration depends weakly on halo mass. 



much more strongly on the orbital anisotropy than on the 
actual density distribution of the parent halo. Since there is 
no obvious reason why the anisotropy of the orbits of sub- 
haloes should depend on halo mass, our assumption that 
the average, instantaneous mass loss rate of substructure 
depends only on m/M should be sufficiently accurate. 



3.1 The average mass loss rate 

During the evolution of the system the parent mass, M, will 
increase due to merging and accretion, whereas the subhalo 
mass, m, will decrease due to the effects discussed above. 
We postulate that in a steady-state halo, for which M = 
AM /At — 0, the instantaneous, fractional mass loss rate of 
a dark matter subhalo is given by rn/m = f(ip) with f(x) 
an arbitrary function (0 < x < 1), to be determined below. 
In what follows we assume, for simplicity, that f(x) is well 
described by a power-law, and write 

m=--V C (1) 

T 

Here r is a characteristic time scale (in Gyr), and £ is an ad- 
ditional free parameter that specifies the mass dependence 
of the subhalo mass loss rate. The negative sign is to empha- 
size that m is expected to decrease with time. For a subhalo 
embedded in a static parent halo (M — 0) this yields 

{ra;exp(— t/r) if C = 

m ,[l + C^ C (t/r)]" 1/C otherwise (2) 

where we have used the boundary condition ip(t = 0) = 
i>i = nii/M. 

We emphasize at this stage that the power-law form of 
the average mass-loss rate has no physical motivation. We 
choose it purely for simplicity. Note that m has to capture 
the effects of both dynamical friction and tidal stripping. 
One might therefore expect that the mass loss rates of indi- 
vidual subhaloes differs significantly from the simple power- 
law form adopted here. However, recall that we use rh to de- 
scribe the average mass loss rate, which is not necessarily of 
the same form as that of individual subhaloes. Furthermore, 
as we show below, it does seem able to naturally explain 
the subhalo statistics found in numerical simulations, with- 
out the need for a more complicated functional form. Nev- 
ertheless, a detailed comparison against the average mass 
loss rates obtained from numerical simulations is required 
to check whether our power-law form is truly appropriate. 

One naturally expects the characteristic time scale r to 
be related to the dynamical time, td yn , of the parent halo. 
As shown in Appendix A, in the idealized case of homolo- 
gous haloes, the mass loss rate of a subhalo on a circular 
orbit can indeed be written in the form (1) with r = tdyn- 
Since td yn oc p~ 1//2 , and since the average density of a dark 
matter halo is a function of redshift, we thus expect that 
r — t(z). The average density of a virialized dark matter 
halo at redshift z is given by p(z) = A v i r (z) p CI i t (z). Here 
pcrit(z) = 3H 2 (z)/8tvG is the critical density for closure, 
and A v i r (,z) is a cosmology-dependent quantity for which 
we use the fitting function of Bryan & Norman (1998). To 
take proper account of the expected redshift dependence of 
the characteristic time scale for subhalo mass loss, we write 
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Figure 3. Left-hand panel: The solid histogram indicates the average, evolved SHMF for a parent halo with Mq = 10 15 /i Mq, obtained 
from 2000 merger trees with to = 0.13 Gyr and £ = 0.36. With these parameters the resulting SHMF best matches those of Gao et 
al. (2004) and De Lucia et al. (2004), shown as dashed and dot-dashed lines, respectively. Note that the model reveals a high-mass cut-off. 
The dotted histogram indicates the unevolved SHMF (i.e., without subhalo mass loss) and is shown for comparison. Right-hand panel: 
Same as left-hand panel, except that here we plot the mass fraction of dark matter subhaloes. Note that d/ s /dln(m) of the evolved 
subhaloes is very flat, indicating that low mass subhaloes contain a significant fraction of the total subhalo mass. 



with to a free parameter that expresses the characteristic 
time scale for subhalo mass loss at z = 0. 



3.2 Evolution of the population of subhaloes 

Although the subhalo mass loss rate in a static parent halo 
is a meaningful concept from a physical point of view, in 
reality parent haloes themselves evolve due to merging and 
accretion. In order to take this into account we utilize the 
discrete time stepping of our merger trees. At the beginning 
of each time step the parent halo is assumed to increase its 
mass through (instantaneous) mergers (M > 0, m = 0), 
while during the period in between two merger events we 
set M — and evolve m(t) according to eq. (2). The exact 
procedure is illustrated graphically in Fig 2: At t — ti halo 
1 (with three subhaloes) and halo 2 (with two subhaloes) 
merge. Since Mi > M2, halo 1 is considered the new parent 
halo, with halo 2 as a subhalo. In addition, the subhaloes 
of Mi are preserved, and are considered subhaloes of the 
new, merged halo. The two subhaloes of 2, however, are no 
longer considered (i.e., we do not follow the evolution of 
sub-subhaloes) . From time ii to t2, which is when the next 
merging or accretion event occurs, the subhaloes evolve ac- 
cording to our mass loss rate, i.e., eq. (2) with t = t% — t%, 
rm = m(ti), and t = 7"(ii). This procedure, hereafter re- 
ferred to as the 'Monte-Carlo method', yields, at each red- 
shift, the evolved SHMF. In addition, we also register for 
each subhalo the time of merging, t m , as well as its mass at 
that time, m(t m ). The abundance of these progenitor haloes 



as function of their mass, m(t m ), is hereafter referred to as 
the unevolved SHMF. 

In order to calibrate our model, we tune the free pa- 
rameters ( and to such that the SHMF of parent haloes 
with Mo = W^h" 1 Mq matches the subhalo mass func- 
tions of Gao et al. (2004) and De Lucia et al. (2004). Al- 
though these two SHMFs, obtained from independent nu- 
merical simulations, are very similar, the agreement is not 
perfect. Including the results of Tormen et al. (2004) we 
estimate the accuracy of the absolute normalization to be 
about 20 percent, and caution the reader that the absolute 
normalization of our results is therefore uncertain by a simi- 
lar amount. Nevertheless, our relative normalizations, which 
are the main topic of interest here, should not be effected 
by this. A more robust absolute normalization will have to 
await a larger sample of high-resolution simulations, and a 
more detailed investigation of numerical resolution effects. 

The dotted histogram in the left-hand panel of Fig. 3 
plots the average unevolved SHMF obtained from 2000 
merger trees for a parent halo of Mo = 10 15 ft _1 Mq, and 
is shown for comparison with the evolved SHMF (solid his- 
togram). The latter is obtained from the same 2000 merger 
trees using the method described above with to = 0.13 Gyr 
and £ = 0.36. These are the parameters for which we ob- 
tain the best-fit to the subhalo mass functions of Gao et 
al. (2004) and De Lucia et al. (2004), shown as dashed and 
dot-dashed curves, respectively. The agreement with these 
SHMFs obtained from numerical simulations is very satis- 
factory, except for a high-mass cut-off in our model, which 
is not accounted for in the simple power-law fits to the pub- 
lished SHMFs of Gao et al. (2004) and De Lucia et al. (2004) . 
Detailed tests have shown that the location of this high- 
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Figure 4. Parent halo mass dependence of the SHMF. Upper left-hand panel: Average SHMF for parent haloes of three different masses, 
as indicated. Less massive parent haloes contain fewer subhaloes at any given mass ratio m/M . For comparison, the dot-dashed histogram 
indicates the unevolved SHMF (which is virtually indistinguishable for parent haloes of different mass). The thin, solid lines are the 
SHMF fitting functions described in Section 7. Upper right-hand panel: Distributions of the ratio of present day mass to mass at time 
of accretion, mo/m(t m ), for all subhaloes in parent haloes with Mo = I0 15 h~ 1 Mq (solid histogram), Mo = 10 13 /i — 1 Mq (dashed 
histogram), and Mo = lO 11 ^ - 1 Mq (dotted histogram). Note that subhaloes in more massive parent haloes have, on average, lost a 
smaller fraction of their mass. Lower left-hand panel: The average subhalo mass fraction, (f s ) (averaged over 2000 merger trees), as 
function of parent halo mass: subhaloes in more massive parent haloes contain a larger mass fraction. Lower right-hand panel: SHMFs 
scaled per unit parent halo mass. Note that with this scaling, the mass dependence is completely removed (except for the high-mass 
cut-off, which simply reflects that m/M < 1). The thin, solid line labelled G04 indicates the results (eq. [6]) obtained by Gao et al. (2004) 
from high-resolution, numerical simulations, and is in excellent agreement with our model predictions. 



mass cut-off is robust to changes in to and/or (. The former 
mainly influences the absolute normalization, while the lat- 
ter controls the slope at small m/M. 

The right-hand panel of Fig. 3 plots 



dlnm dlnm 

which indicates the mass fraction of the parent halo that is 
associated with subhaloes of mass m. Most of the present 
day halo mass originates from progenitors (i.e., unevolved 



subhaloes) with m(t m ) ~ OTA/o. In fact, if one ignores 
all progenitors with m(t m )/M~o < 1CF 4 one misses only a 
negligible fraction of the entire mass. This, however, is not 
the case for the evolved SHMF. Here d/ s /dlnm is remark- 
ably flat, indicating that even the evolved subhaloes with 
m/M < 10~ 4 contribute a significant fraction of the to- 
tal subhalo mass. Therefore, it is important to always in- 
dicate the range of m/M considered when quoting subhalo 
mass fractions, and care is required when comparing sub- 
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halo mass fractions from different simulations with different 
resolutions. 



4 THE MASS DEPENDENCE OF THE 
SUBHALO MASS FUNCTION 

Having calibrated our mass-loss rate by matching the sub- 
halo mass function obtained from numerical simulations 
for haloes with Mo = 10 15 /i -1 M©, we now investigate 
what the Monte-Carlo method predicts for different par- 
ent halo masses. The upper, left-hand panel of Fig. 4 plots 
the SHMFs obtained using the Monte-Carlo method with 
to = 0.f3 Gyr and £ = 0.36 for three different parent halo 
masses, as indicated (each SHMF is averaged over 2000 
merger trees). For comparison, we also show the unevolved 
subhalo mass function, which is virtually identical for all 
three halo masses. The evolved SHMFs, however, are clearly 
mass-dependent with a normalization that decreases sys- 
tematically with decreasing halo mass. These findings are 
in good agreement with those of Gao et al. (2004), based 
on numerical simulation that span three orders of magni- 
tude in parent halo mass, and strongly argue against earlier 
claims for a universal subhalo mass function (e.g., Moore et 
al. f998; De Lucia et al. 2004). 

The mass dependence of the normalization of the 
evolved SHMF is simply a reflection of the fact that less 
massive haloes form earlier, thus providing more time for 
mass loss to operate. This is illustrated in the upper, 
right-hand panel of Fig 4, which plots the distributions of 
mo/m(tm) for parent haloes of Mo = lO 11 ^ 1 M© (dot- 
ted histogram), Mo = 10 13 /i _1 M© (dashed histogram), 
and Mo = 10 1B /i _1 M© (solid histogram). These distribu- 
tions, clearly show that subhaloes in less massive parent 
haloes have, on average, lost a relatively larger fraction of 
their mass since they were accreted. The distributions of 
mo/m(t m ), however, are very broad; while some subhaloes 
have only lost a negligible fraction of their initial mass (ei- 
ther because they were accreted relatively late, or because 
they had small, relative masses to begin with), others have 
lost more than 99.9 percent of their mass since their time of 
accretion. 

4.1 The mass fraction in subhaloes 

To quantify the mass dependence of the SHMF we con- 
sider the mass fraction of dark matter subhaloes with m > 
10- 4 M: 

f.= f V^ch/> (5) 

The lower limit of 10~ 4 reflects the effective mass resolu- 
tion of our merger trees (see Section 3). Throughout this 
paper all subhalo mass fractions therefore only take account 
of subhaloes with masses above this resolution limit. As dis- 
cussed above, the total subhalo mass fraction can easily be 
a factor two larger than this. 

The lower, left-hand panel of Fig. 4 plots (f s ) at z = 
(where the average is taken over 2000 merger trees) as func- 
tion of parent halo mass. This nicely illustrates the mass 
dependence of the SHMF indicated above, namely a system- 
atic increase of (f s ) with increasing parent halo mass. From 




<5,=f s /<f s > " 1 

Figure 5. Distributions of Sf (eq. [7]), obtained from 2000 in- 
dependent merger trees, for three different parent halo masses as 
indicated. Note the strong skewness, and the relatively large dis- 
persion. The vertical line indicates the average for which Sf = 0. 
See text for detailed discussion. 

the scale of galaxy sized haloes (Mo — 10 12 /i _1 M©) to that 
of massive clusters (M ~ lO 15 ^ 1 M©) we find that (f B ) 
increases by about a factor three, in reasonable agreement 
with Gao et al. (2004) and Kang et al. (2004). Note that 
(f s )(M) seems to asymptote to a non-zero value of ~ O.Of 
for low mass parent haloes, which is simply a reflection of 
the finite age of the Universe; only in the limit where the 
formation time of a halo is infinitely long ago will there be 
sufficient time to wipe out all substructure. 

4.2 Subhalo mass functions per unit halo mass 

The lower, right-hand panel of Fig. 4 shows the subhalo mass 
functions for six different parent halo masses (each obtained 
using ro = 0.f3 Gyr and £ = 0.36, and averaged over 2000 
merger trees), but this time normalized in a different way. 
Following Gao et al. (2004) we divide the total number of 
subhaloes in each bin by the total mass of all the parent 
haloes (in units of 10 10 /i -1 M©) to obtain the subhalo abun- 
dance per unit parent halo mass. These abundances are plot- 
ted as function of the actual subhalo mass, m, rather than 
the scaled mass, m/M. With this particular normalization, 
the subhalo mass functions of different parent halo masses 
agree extremely well (except for the high-mass cut-off). The 
thin, straight line corresponds to 

*L = io-" (-^—Y'\m^ (6) 

dm yh- 1 M© J w 

which is the best-fit subhalo abundance per unit halo 
mass (ignoring the high mass cut-off) obtained by Gao et 
al. (2004). As can be seen, our results are in excellent agree- 
ment with those of Gao et al. (2004), lending strong support 
for our simple model. 
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Figure 6. Distributions of parent halo formation times, tf orm , (upper left-hand panel) and of the mass fraction AM(1)/Mq accreted 
in the last 1 Gyr (lower left-hand panel). Results are shown for three different parent halo masses, as indicated. The right-hand panels 
show, for a parent halo mass of Mo = 10 ls h~ 1 Mq, how tf orm and AM(1)/Mq correlate with Sf. The relatively tight relation between 
Sf and AM(l)/Mo indicates that the subhalo mass fraction of individual haloes depends mainly on their accretion history in the last 
~ 1 Gyr. 



5 SCATTER IN SUBHALO MASS FUNCTIONS 

Thus far we only focussed on the average SHMFs, where 
the average is taken over all orbital configurations and over 
many mass accretion histories (hereafter MAHs). However, 
since there is considerable scatter in MAHs of parent haloes 
of the same mass, and since individual haloes may have sig- 
nificantly different orbital distributions for their subhaloes, 
one expects a relatively large halo-to-halo variation in the 
SHMF. Here we use the Monte-Carlo method to obtain an 
estimate of this scatter. Since this method implicitly aver- 
ages over all orbital configurations, we can only address the 
halo-to-halo scatter due to variance in the MAHs. Our esti- 
mates of the amount of scatter are therefore to be considered 
lower-limits. 



Fig. 5 plots the distributions of 

_ fs - (fs) 
- (fs) 



(7) 



obtained from 2000 independent MAHs (merger trees). Note 
that these distributions are extremely skewed (not surpris- 
ing, given that f s > 0), and fairly broad. This indicates that 
the SHMF obtained from a small number of haloes, as is 
typically the case with current simulations, may not be an 
accurate representation of the true, average mass function. 
This explains, at least partially, why it is so difficult to in- 
fer from numerical simulations whether or not the SHMF 
depends on parent halo mass; only when averaged over a 
sufficiently large number of parent haloes will such a trend 
become evident. 



© 2000 RAS, MNRAS 000, 1-14 



The Mass Function and Average Mass-Loss Rate of Dark Matter Subhaloes 9 




1.5 



A i 



0.5 - 



log[m/M]^-4 

log[m/M]^-3 
i I i i i l_ 



10 12 14 

log[M] (h"«M ) 




12 14 
log[M] (h"«M ) 



Figure 7. The upper panels plot the distributions P(N), where N is the number of subhaloes with log[i/>o] > — 4 (upper left panel) 
and log[«/>o] > — 3 (upper right panel), respectively. Results, obtained from 2000 merger trees, are shown for three different parent halo 
masses, as indicated. The lower panels plot (N(N - l)) 1 / 2 /(A r ) (lower left panel) and (N(N - 1)(N - 2)) 1 / 3 /(N) (lower right panel), 
which express the second and third moments of P(N) in the units of that of a Poisson distribution with the same (TV). See text for a 
detailed discussion. 



The upper, left-hand panel of Fig. 6 plots the distri- 
butions of the parent halo formation times, tform, defined 
as the lookback time at the redshift where the mass of 
the most massive progenitor reaches half the present day 
mass. Although these distributions are very broad, there is 
a clear mass-dependent trend in that more massive haloes 
form later. As discussed in Section 4, this explains why (/ s ) 
increases with parent halo mass; in systems that form later, 
the subhaloes have less time to experience mass loss. It 
therefore seems natural that the scatter in iform is the direct 
source of scatter in f 3 . However, as evident from the upper, 
right-hand panel of Fig. 6, the correlation between tf or m and 
5f is, in fact, surprisingly weak. 

We find a must stronger correlation between Sf and 
AM(1)/M (lower right-hand panel of Fig. 6), with AM(1) 



the mass that has been accreted by the parent halo in the 
last 1 Gyr. This is easy to understand. Since the charac- 
teristic time scale for subhalo mass loss is relatively short 
(tq = 0.13 Gyr) compared to the typical halo formation 
time, the subhalo mass fraction in individual systems is dom- 
inated by the mass that was accreted relatively recently; we 
find that the scatter between <5/ and AM(t)/Mo is mini- 
mized for t ~ 1.0 Gyr, which is the value adopted here. 
A similar conclusion was reached by Gao et al. (2004) who 
found that subhaloes are typically recent additions to their 
parent haloes, substantially more recent, in fact, than typical 
dark matter particles (see also ZB03). Thus, whereas the av- 
erage SHMF depends strongly on formation time, the SHMF 
of an individual halo simply reflects its accretion history in 
the last ~ 1 Gyr. Although this may seem contradictory, 
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shown as a dot-dashed histogram. Note that haloes at high redshift contain, on average, more subhaloes than haloes of of the same mass 
at lower redshifts. As in Fig. 4 the thin, solid lines are the SHMF fitting functions described in Section 7. Right-hand panel: Redshift 
dependence of the average subhalo mass fraction, (f s ), for parent haloes of 5 different masses. Each curve is labelled by M(z) (in h~ 1 'Mo). 



it is easy to understand when looking at the distributions 
of AM(1)/Mq, which are shown in the lower-left panel of 
Fig 6 for parent haloes of three different masses. As is evi- 
dent, more massive haloes have, on average, accreted more 
mass recently, which reflects their relatively later formation 
times, and which is responsible for the mass-dependence of 
the average SHMF. 

In addition to the scatter in the subhalo mass fraction 
f s , we also investigate the scatter in the number of subhaloes. 
The upper two panels of Fig. 7 plot the distributions P(N), 
where N is the number of subhaloes with log [ipo ] > —4 
(upper left panel) and log^o] > ~3 (upper right panel), 
respectively. Results are shown for three different parent 
halo masses, as indicated. As above, these probability dis- 
tributions reflect 2000 independent MAHs (per parent halo 
mass). These two panels, once again, clearly demonstrate 
that the subhalo mass function is not universal, but instead 
depends strongly on parent halo mass: more massive haloes 
contain more subhaloes above a give ^-threshold. The lower 
two panels plot M 2 = (N(N - l)) 1/2 /(N), related to the 
second moment, and M 3 = {N(N - 1)(N - 2)) 1/3 /{N), re- 
lated to the third moment, of these distributions, as function 
of parent halo mass M. Note that for a Poissonian distribu- 
tion M2 = M-3 = 1, while distributions that are narrower 
(sub-Poissonian) or broader (super-Poissonian) have Ai < 1 
and M > 1, respectively. Clearly, when considering all sub- 
haloes with log[^o] > — 4 (solid lines), the P(N) are very 
close to Poissonian. However, when only counting the more 
massive subhaloes, with logf^o] > —3, the distributions are 
super-Poissonian, with a clear trend of increasing Ai with 
decreasing parent halo mass. These results are inconsistent 
with Kravtsov et al. (2004), who finds that the number of 
subhaloes in numerical simulations follow Poisson statistics. 
This may reflect a generic problem of the EPS formalism 



used here to construct the merger trees: as shown by various 
authors (Lacey & Cole 1993; Somerville et al. 2000; van den 
Bosch 2002; Wechsler et al. 2002), the halo formation times 
predicted by EPS are systematically offset from those ob- 
tained from numerical simulations. In particular, Somerville 
et al. (2000) found the average mass of the largest progenitor 
to be larger with the EPS formalism than in the simulations. 
This may explain why we find the non-Poissonian nature of 
P(N) to be more pronounced for more massive subhaloes. 
Although the merger trees extracted from numerical simu- 
lations have their own problems, we caution that the scatter 
issues discussed here are probably less reliable than the av- 
erage mass trends. 



6 REDSHIFT EVOLUTION OF THE SUBHALO 
MASS FUNCTION 

The left-hand panel of Fig. 8 shows the average subhalo 
mass functions for parent haloes of the same mass, M(z) = 
10 15 /i _1 Mq, but at different redshifts. At higher redshifts 
parent haloes of the same mass have a larger abundance of 
subhaloes than their counterparts at lower redshifts. This 
is quantified more clearly in the right-hand panel of Fig. 8, 
which shows the redshift dependence of the average subhalo 
mass fraction. The various curves are labelled by the parent 
halo mass M(z) (in h~ x Mq). In all cases {f B )(z) increases 
with redshift, though with a rate, d(f s )/dz, that decreases 
monotonically. Roughly speaking, the subhalo mass fraction 
at z = 1 is about twice as large as that of a halo with the 
same mass at z = 0. 

The subhalo mass fraction of a given halo at redshift 
z is a trade-off between the time scale, f acc , on which new 
subhaloes are being 'accreted' by the parent halo, and the 
time scale, r, of subhalo mass loss. The latter evolves with 
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Figure 9. The dependence of (/ s ) (top panel), a (middle panel), 
and /3 (bottom panel) on the mass ratio M/M*. Solid and open 
circles correspond to parent haloes in a ACDM concordance cos- 
mology and a SCDM cosmology, respectively. Both the normal- 
ization and the shape parameters a and j3 are tightly correlated 
with M/M* , indicating that the latter is the main parameter de- 
termining the average SHMF of dark matter haloes. The solid 
lines are the best-fit functions (9), (13), and f) = 0.13 discussed 
in the text. 



redshift as described by eq. (3), and therefore was shorter 
in the past. The former depends on the detailed MAH, and 
is thus a function of both redshift and parent halo mass. 
In the limit where t &cc <C t, subhalo mass loss is negligible 
and f s will increase with time. The opposite limit, in which 
iacc 3> t, is equivalent to that of subhalo mass loss in a 
static parent halo. In this case, f s will decrease with time. 
Since the subhalo mass fraction always decreases with time, 
the time scale for subhalo mass loss is always smaller than 
that of mass accretion; r < t acc . 



7 AN ANALYTICAL FITTING FUNCTION 
FOR THE SUBHALO MASS FUNCTION 

Since the unevolved SHMF is virtually independent of M, 
z, or even of the cosmological parameters (see Lacey & Cole 
1993 and ZB03 for detailed discussions), and since we have 
adopted a universal, average subhalo mass loss rate, the av- 
erage, evolved SHMF simply depends on the halo formation 
time. This suggests that the mass, redshift, and cosmology 
dependence of the average SHMF can be written as a sim- 
ple one-parameter dependency on the mass ratio M/M*, 
with the characteristic non-linear mass M*(z) defined by 
a(M* , z) = 5 c (z). Here a 2 (M, z) is the mass variance of the 
smoothed density field at redshift z and 



S c (z) = 0.15(T2tt) 2/3 [n m (z)}° 



(8) 

Navarro 



is the critical threshold for spherical collapse (e.g.. 
et al. 1997). 

The top panel of Fig. 9 plots the average subhalo 
as function of M/M* (solid circles), 
for six different parent halo masses 
10, 11, 15) at six different redshifts 



mass fraction, (/ s ), 
Results are shown 
(loglM/h- 1 M c . 



iQj = 

(z = 0,1, ..5). Although this yields values of M/M* up 
to 10 s , we caution that systems with M/M* •> 10 4 are ex- 
tremely rare. The open circles indicate the results for the 
same halo masses and redshifts, but obtained for a SCDM 
cosmology with Q m = 1.0, Qa = 0.0, h = 0.5, and cr& = 0.7. 
All these different haloes follow a tight relation between (/„) 
and M/M*, which is well fitted by 



log[(/„>] = V0-4(log[Af/Af*] + 5) - 2.74 



(9) 



indicated by the solid line. This indicates that, as expected, 
the average subhalo mass function is completely specified 
by the mass ratio M/M* . To quantify this further we fit the 
average subhalo mass functions with a Schechter function of 
the form 



7 



dn 



dlm/> /3T(l-a) 



exp 



(10) 



(cf. Vale & Ostriker 2004). Here (3M is a characteristic mass, 
such that for m 3> /3M the SHMF reveals an exponential 
decline, and 7 is the total subhalo mass fraction, i.e., 



7 = 



(11) 



Some of these fits are shown as thin, solid lines in Fig. 4 
(upper left-hand panel) and Fig. 8 (left-hand panel), and 
in general match the model SHMFs extremely well. Note 
that in practice we don't fit the actual SHMF, dn/dlni/), but 
rather the corresponding dfs/dimp, treating a and f3 as free 
parameters. The normalization, 7, is not treated as a free 
parameter, but is fixed by requiring to match the subhalo 
mass fraction f s , which implies 



7 = 



P(l-a,l//?)-P(l-a,10-7/3) 



(12) 



with P(a, x) the incomplete Gamma function. 

The middle and bottom panels of Fig. 9 plot the best- 
fit a and j3 for each of the 72 SHMFs (2 cosmologies, 6 
masses, 6 redshifts) as function of the mass ratio M/M* . As 
with the average subhalo mass fraction, a and (3 are both 
tightly correlated with the parent halo mass in units of the 
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characteristic non-linear mass; the power-law slope a scales 
roughly linearly with log(M/M*), which is best fit by 

a = 0.966 - 0.028 log(M/M*), (13) 

while the parameter (3 is best fit by f3 — 0.13. 

We thus have obtained an extremely simple recipe to 
compute the average subhalo mass function for a parent halo 
of any mass, at any redshift, and for any cosmology: compute 
the characteristic non-linear mass M*(z), and use eq. (9) 
and (13) to obtain both f s and a. The normalization 7 then 
follows from (12), which, together with f3 = 0.13, completely 
specifies the SHMF. 



8 CONCLUSIONS 

We combined merger trees of dark matter haloes, con- 
structed using the EPS formalism, with a simple prescription 
of the average mass loss rate of dark matter subhaloes, to 
compute subhalo mass functions. We calibrated the subhalo 
mass loss rate by matching the SHMFs of massive haloes 
with Mo = 10 15 /i _1 Mo obtained from high-resolution, nu- 
merical simulations. Under the assumption that this average 
mass-loss rate only depends on redshift and on the mass ra- 
tio of sub- and parent halo, m/M, this method allows us 
to make detailed predictions for the mass and redshift de- 
pendence of the SHMF, and to investigate the halo-to-halo 
variance. 

Our main conclusions are: 

• Contrary to previous claims, the subhalo mass function 
is not universal. Instead, both the slope and the normaliza- 
tion depend on the ratio of parent halo mass, M, to charac- 
teristic non-linear mass, M* . This simply reflects a halo for- 
mation time dependence, in which parent haloes that form 
earlier have a lower subhalo mass fraction, because there is 
relatively more time for subhalo mass loss to operate. 

• When the subhalo mass function is normalized by the 
mass of the parent halo, the abundance of subhaloes is uni- 
versal, in excellent agreement with the numerical simulations 
of Gao et al. (2004) and Kravtsov et al. (2004). 

• The subhalo mass function of an individual halo de- 
pends most strongly on its accretion history in the last ~ 1 
Gyr. This indicates, as previously shown by ZB03 and Gao 
et al. (2004), that the population of the more massive dark 
matter subhaloes is, at any time, relatively young. 

• The dependence of SHMF on the recent accretion his- 
tory introduces a large halo-to-halo scatter in the SHMFs, 
with a distribution of subhalo mass fractions that is strongly 
skewed towards large values. 

• The average subhalo mass function of dark matter 
haloes is a one-parameter family, depending only on the 
mass ratio M/M*. We have provided simply fitting func- 
tions that allow one to compute the average SHMF for a 
parent halo of any mass, at any redshift, and for any cos- 
mology. 

While this paper was being refereed, a paper appeared 
by Zentner et al. (2004) which uses a similar semi-analytical 
model as in ZB03 and TB04 to compute subhalo statistics. 
Using a proper integration of individual orbits, and taking 
detailed account of dynamical friction and tidal stripping, 



these authors reach conclusions that are in excellent agree- 
ment with the simplified method presented here. In particu- 
lar, they find that (i) the subhalo mass function is not uni- 
versal but scales with halo mass, (ii) the subhalo mass func- 
tion is most sensitive to the most recent accretion history 
of the parent halo, and (iii) the distribution of the number 
of subhaloes per parent halo is super-Poissonian. The good 
agreement of this more sophisticated model with that pre- 
sented here, provides further support for our 'orbit averaged' 
approach. 

These results have a number of important, astrophys- 
ical implications. For example, the prediction that the av- 
erage subhalo mass fraction in galaxy sized haloes is a fac- 
tor three lower than in cluster-sized haloes has important 
implications for the magnitude of the claimed substructure 
crisis (Moore et al. 1999; Klypin et al. 1999b: D'Onghia & 
Lake 2004), for the flux-ratio statistics of multiply lensed 
quasars (e.g., Chiba 2002; Metcalf & Madau 2001; Bradac 
et al. 2002; Dalai & Kochanek 2002), and for the build-up of 
the galactic halo (Helmi, White & Springel 2003). Further- 
more, the redshift dependence of the subhalo mass fraction 
impacts on the survival probability of fragile structures in 
dark matter haloes, such as tidal streams and/or galactic 
disks (Toth & Ostriker 1992; Taylor & Babul 1991). The 
subhalo mass functions derived here may also be used in 
combination with the so-called halo model (see Cooray & 
Sheth 2002 and references therein) to give a full statistical 
description of the distribution of dark matter haloes down 
to the level of subhaloes. This will proof especially fruitful 
in combination with the conditional luminosity function for- 
malism developed by Yang et al. (2003) and van den Bosch 
et al. (2003), allowing for a detailed, statistical description of 
the relation between light and mass (see also Vale & Ostriker 
2004). Finally, the average mass loss rates derived here may 
proof useful in semi-analytical models for galaxy formation, 
where a proper treatment of the evolution of subhaloes is ex- 
tremely important (Springel et al. 2002; Benson et al. 2002; 
Kang et al. 2004). 

Finally we point out that, although the Monte-Carlo 
method presented here is nice and simple, it is important 
to be aware of its potential shortcomings. For example, the 
accuracy of the absolute normalization of our model is only 
as good as that of the SHMFs used for its calibration, which 
we estimate to be about ~ 20 percent. In the numerical sim- 
ulations used for this calibration, the masses of the parent 
haloes are defined as the masses within a sphere of density 
200 times the critical density at redshift zero. Therefore, we 
have implicitly assumed that the masses in the EPS for- 
malism used to construct our merger trees are defined in 
the same way. Since it is still unclear what the proper in- 
terpretation of these EPS masses is (see White 2002 for a 
detailed discussion), this 'definition' may introduce an ad- 
ditional uncertainty in the absolute normalization of our re- 
sults. The fact that the method to construct merger trees is 
not without its own shortcomings (e.g., SK99; TB04; Sheth 
& Tormen 1999; Benson, Kamionkowski & Hassani 2004) 
may have additional implications for the accuracy of our re- 
sults. Furthermore, we have ignored the weak dependence of 
halo concentration on halo mass, which may cause a (weak) 
dependence of the average subhalo mass loss rate on M in 
addition to tp (see e.g., ZB03). Finally, we have ignored any 
possible effect due to subhalo-subhalo mergers. Although the 
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excellent agreement between our results and those of Gao et 
al. (2004) suggest that none of these effects have a strong 
impact, large, high-resolution numerical simulations are re- 
quired to further test both the validity of our approach as 
well as the accuracy of our results. 
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APPENDIX A: THE MASS LOSS RATE OF 
SUBHALOES ON CIRCULAR ORBITS 

Consider a subhalo with density distribution p s (r) on a cir- 
cular orbit in a parent halo with density distribution p P (r). 
For simplicity, we assume that both p s (r) and p P (r) are sin- 
gular, isothermal spheres. 

In the absence of tidal heating, the mass loss rate of the 
subhalo is given by 



dm _ dm dr«d 
dt dr t id dt 



(Al) 



with r t id the instantaneous tidal radius of the subhalo. We 
define this tidal radius as the radius of the subhalo where 
its density is equal to that of the parent halo at the orbital 
radius r or b of the satellite: p s (r t id) = p P (r or b). Since p s (r) 
and p p (r) are scale- free we have that 



1 dr tid 



1 dr orb 



rtid dt r orb dt 



(A2) 



The evolution of r orb is governed by dynamical friction, and 
is given by 



dr, 



orb 



dt 



-0.428 #^-lnA 

Vc r OT b 



(A3) 
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(Binney & Tremaine 1987) with V c the circular velocity ol 
the parent halo and InA the Coulomb logarithm, which, to 
leading order, is just a lunction of the mass ratio m/M (Bin- 
ney & Tremaine 1987) . 

Using that dm/dr t id = 47rp s (r t id) r 2 id and V c 2 = 
GM/r v i r , with r vlI the virial radius of the parent halo, we 
obtain that 



_«-p.(r tld )Vc-^ 5f InA (A4) 



, -—( I™.) HL InA (A5) 
\r r h ) M K ' 



dm , , T T rl 

' orb 

Defining the dynamical time as £dyn oc r v ; r / V c and using that 
m tx p s (r t id)rl id , the subhalo mass loss rate can be written 

as 

dm ^ m ( r v ; r \ 2 m 

\ ' orb ' 

Therefore, when averaging over all possible circular orbits, 
i.e., all possible ratios r v ir/V rb, one obtains an average mass 
loss rate for which 

dm m , , , „. 

— cx-— g (m/M) (A6) 

at tdyn 

with g(x) an arbitrary function. This is the basic form for 
the average mass loss rate adopted in this paper. Note that 
since tdyn oc p -1 / 2 , and since the average density of dark 
matter haloes change with redshift, eq. (A6) automatically 
implies a redshift dependence (see Section 3.1). 



